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Abstract. The analysis of diffusive energy spreading in quantized ehaotic driven 
systems, leads to a universal paradigm for the emergence of a quantum anomaly. 
In the classical approximation a driven chaotic system exhibits stochastic-like 
diffusion in energy space v^ith a coefficient D that is proportional to the intensity 
£^ of the driving. In the corresponding quantized problem the coherent transitions 
are characterized by a generalized Wigner time , and a self-generated (intrinsic) 
dephasing process leads to non- linear dependence of _D on e^. 

A major theme in mechanics concerns the response of a system to a driving source 
/(i), given that the interaction term is Hi„t = f{t)V. This leads to the well known 
framework of linear response theory (LRT) with its celebrated fluctuation-dissipation 
relation. Below we assume a stationary driving source which is characterized by a 
power spectrum S{uj) — FT[(/(i)/(0))], where FT stands for Fourier transform. In 
the absence of driving the stationary fluctuations of the system are characterized by 
the spectral function C{uj) = FT[{V{t)V{0))]. In the presence of driving the main 
three effects are: the decay of the initial preparation; the spreading and eventually the 
fusion in energy space; and the associated heating. 



1. LRT and Kubo 

Strict LRT behavior means that the diffusion in energy space [I] and the related 
absorption coefficient [21 131 H] are linear functional of the spectral function S{io). 
Specifically, the Kubo formula for the diffusion coefficient in energy space is 

^ r°° > > 

D = - u'^duj C{uj)S{uj) (1) 



It follows that the diffusion is proportional to the intensity of the driving as 
defined below. We are going to consider on equal footing driving by a quasi-constant 
perturbation f{t) ^const, and quasi-linear DC driving with f(t) ^const. The notation 
"~const" means that it is constant over large time intervals of duration , with some 
characteristic RMS value that we call e. Accordingly the associated spectral functions 
is 

S{uj) = e^uj-^S^iuj) (2) 

where S.y{uj) = {■^/'k)/{uP +7^) with 7 = 1/t^p, while the spectral exponent is ct=0 for 
quasi-constant perturbation, and <j~2 for quasi linear DC driving. 
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2. Paradigms 

We are looking for circumstances where the Kubo formula is not applicable. This 
means that either D is still proportional to the strength of the driving but with an 
anomalous coefficient, or more generally D might depend in a non-linear way on the 
strength of the driving. Paradigms for non-LRT response with which we are familiar 
are: (i) Classical non-LRT route that follows from the Kolmogorov-Arnold-Moser 
scenario as in the driven (kicked) rotator problem [5]; (ii) Quantum semi-LRT due 
to sparsity and textures that characterize the perturbation matrix ; (iii) Quantum 
corrections due to dynamical localization effect [71 [5] . (iv) Quantum absorption due 
to Landau-Zener transitions between neighboring energy levels [3; (v) Quantum non- 
perturbative anomalies that are associated with having a finite spectral bandwidth ^ . 

3. Scope and main observation 

One should notice that all the mentioned non-LRT paradigms above become irrelevant 
if we consider the continuum limit of a universal quantized chaotic system. By 
definition "continuum" means that the Heisenberg time (see definition later) can be 
taken as infinite, while "universality" means that the correlation time (see definition 
later) can be taken as zero. It should be clear that our considerations apply to strictly 
chaotic systems that do not have mixed phase-space dynamics. In the present work 
we show that even with all these assumptions and exclusions there is still room for 
a novel manifestation of quantum mechanical anomalies in the response to external 
driving. A central role is played by the generalized Wigner time which characterizes 
a coherent spreading process in energy space, and by what we call intrinsic dephasing 
time tl^"^ . 

Beyond any technical details it is important to realize that for a universal 
chaotic system, as assumed in Random Matrix Theory (RMT) studies, the Kubo 
formula of LRT can be deduced merely via dimensional analysis. The non-Ohmic 
generalization of this statement [Eq. ([TOl) ]. as established in this communication, 
implies a universal quantum anomaly in the response characteristics of non-Ohmic 
systems. This prediction has potential applications e.g. with regard to the rate of 
heating of cold atoms in vibrating traps. 

4. Modeling 

In a system that is described by a time dependent Hamiltonian T-L[R] with R — 
Ro+f{t) the transitions between the adiabatic energy levels En are induced by the 
perturbation matrix Vnm = id'H/dR)nm- The spectral function that characterizes the 
fluctuations of V in the absence of driving is 

C{u) = \Vn,J'2nS L - ^^^^) (3) 

n ^ ^ 

with implicit average over no as determined by the energy window of interest. We 
assume below that 

(7(0;) = 27r|w|''«-i iov luq < \uj\ < co,, (4) 

and distinguish between the Ohmic (50=1)7 subOhmic (0<so<l)j and superOhmic 
(l<,so<2) cases. Without loss of generality, by appropriate rescaling of f{t), we set the 
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prefactor in Eq.Q as 27r. The infrared cutoff ljq — (hg)^^ is the mean level spacing, 
as determined by the density of states. The ultraviolet cutoff us^i is determined by 
the classical dynamics and is known as the bandwidth or as the ballistic version of 
the Thouless energy. The associated time scales are the Heisenberg time = 2TThg 
and the classical correlation time t^, = 2tt/uj^,i. In a later paragraph we define the 
generalized Wigner time te that depends on the strength of the driving. This time 
scale characterizes the coherent spreading process. We assume below mesoscopic 
circumstances for which 



where is the correlation time of the driving source as define with relation to Eq.Q. 
Our interest is in results that remain well defined for i^i — >■ (universal limit) and 
tn oo (continuum limit). The existence of a universal limit is the underlying 
postulate in RMT modeling. We distinguish in the analysis between weak {te > t^) 
and strong [t^ <t^) driving. 

5. DC driving 

The common interest is in Ohmic systems (sq — 1) with quasi-linear DC driving 
(ct = 2), for which Kubo formula gives D = -ne^ . This result is not sensitive to t^p, and 
is independent of the infrared and ultraviolet cutoffs. Our purpose is to generalize 
this result for the case of non-Ohmic fluctuations. We shall see that this requires to 
go beyond LRT. 

The key observation is that the problem of quasi-linear DC driving reduces, with 
some reservations, to the analysis of quasi-constant perturbation. This is done by 
transforming the Hamiltonian into the adiabatic basis where it takes the form 



If we ignore the implicit time dependence of the adiabatic energies and matrix 
elements, then this Hamiltonian is the same as that of quasi-constant perturbation 
but with effective exponent s — ~ 2. In particular quasi- linear driving of an Ohmic 
system corresponds to s = — 1. 

At this point one wonders what is the effect of the residual implicit time 
dependence of the Hamiltonian Eq.([6]). Obviously we should not be too worried about 
the wiggles of the levels, because they take place on a very small energy scale and 
would be of relevance only if we were considering times of the order of the Heisenberg 
time. On the other hand, the variation of the matrix elements cannot be ignored, and 
we shall come back to this issue later. 

6. The generalized Wigner time 

Universality (irrelevance of a;^) is a common built-in assumption in numerous 
"quantum chaos" studies that utilize the standard random matrix ensembles. 
Furthermore a quasi-continuum assumption (irrelevance of loq) is implicit in the 
standard derivations of LRT. If we believe that in the continuum limit (wq — >■ 0) 
there exists a universal limit (a;„i oo) that leads to a generalized response theory, 
then disregarding t^^ the only relevant time scale that might emerge in the dynamics 
is implied by dimensional analysis: 



H 



mesoscopics 



(5) 




(6) 



t, ={hlef/(^-'^ 



(7) 
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One immediately realizes that it is the generalized Wigner time of Ref.|10). For a 
quasi-constant perturbation of Ohmic system it is literally the Wigner time 

t, = (h/r^), for s = 1 (8) 

where Fe = (27r//i)e^ is the Fermi-golden-rule rate of transitions. For a quasi-linear 
driving of Ohmic system it is the breaktime scale that has been introduced in Ref. [4] , 

t, ^ {h^/Df/^, for s = -1 (9) 

where D = 7re^ in derived from the Kubo formula. 

More generally, in the non-Ohmic case, we can associate with the generalized 
Wigner time an energy scale h/te and a diffusion coefficient 

D, = = h-^i^e^, [s = sa-(j] (10) 

At this point one should wonder whether this expression might emerge from the 
analysis of the spreading process in some universal limit. Note that it is only in 
the Ohmic (Kubo) case that De becomes h independent. 

7. Coherent spreading 

If the perturbation matrix in Eq.®, call it Wnm, were strictly time independent, then 
the induced wavepacket dynamics would lead to a steady state, with a saturation 
profile that reflects the local density of states. Specifically, let us assume that the 
system is prepared in the unperturbed state n for which the unperturbed energy En 
is well defined, then in the perturbed basis it has the energy distribution 

P^{E)=J2\{^H?HE-E,) (11) 

We have argued in Ref. [10], following previous studies, that this energy distribution 
has a semicircle-like core that extends within \E—En\ < h/t^, coexisting with outer 
perturbative tails that are determined by the first order expression \ Wnm\'^l{En—Emf' 
for the overlaps. The associated variance is Ai?(oo)^ = a;^,e^ for s > 0, and 
AE{oof = (?iAJ^ for s < 0. 

In the time dependent analysis the steady state profile of Ft (E) is achieved only 
after tg, but the crossover is not necessarily observed in the spreading AE{t), which 



is a second moment calculation. Specifically we get 

AE{t)=suj"J^ for [s>0],[i>y (12) 

AE{t) = e <l''l/2 for [s<Q\\K,<t<Q (13) 

AE{t) = h/t^ for [s<Q\\t>t^] (14) 



For s > the spreading saturates as soon as t > t^i, and to detect the crossover at 
one should look on the survival probability or on percentiles of the distributions as 
described in Ref. 1 10). But for s < the second moment of the evolving distribution 
exhibits the crossover to saturation at t^, and not at t^i. This can be deduced using 
the following simple reasoning. First order perturbation theory implies that the tail 
grows like \Wnmt\'^ within the shrinking interval h'y{t) < \E ~ En\ < h/t. In the outer 
\E—En\ > h/t region the tail is saturated due to recurrences. The lower cutoff 7(t) 
is determined by a self-consistency condition, saying that the integral over the 1st 
order tail of Pt{E) from 7(t) to infinity should be 0(1). This leads to the estimate 
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7(i) — + (£t/?i)^^] Steady state is achieved at when ^{t) ^ l/t. Up to 

this time the second moment of the evolving distribution is dominated by the growing 
piece of the tail leading to Eg. dT^ . This is a diffusive-like growth {AE{t) oc t^/"^) in 
the case of a quasi-linear DC driving. 



8. Diffusion 



The dephasing time tip indicates the crossover from coherent to stochastic spreading 
behaviour. The central limit theorem implies that the long time spreading is diffusive 
with coefficient 

(15) 

where AE{t) is the time dependent coherent spreading that we have discussed in the 
previous paragraph. This result is classical in nature (no h) and it is easily checked 
that it agrees with the Kubo formula Eq.([T]) provided we use Eg. dT^ or Ea. (IT5)) . 
leading to 

D =e^ cj^j t-^ for [s>0] (16) 

D =e^t\^^-^ for [s<0],[<^<t,] (17) 

Both results are oc e^. In particular Eq. (|16|) applies to quasi-constant perturbation 
and is merely the well known hopping estimate for the noise- induced diffusion is system 
with 'localization'. 

Consider a general system with quasi-linear DC driving, such that s < 0. If the 
driving is weak (t^ < te) it is justified to substitute Eq.(fT3|) in Eq. (fT5| . thus getting 
the LRT result Eq. (fT7|) . In particular we note that in the standard case of DC-driven 
Ohmic system (s — —1) we get a <^ independent result. Otherwise the result is 
dependent. For |s| 1 one observes that in the limit ^ oo the LRT result is either 
zero or infinity. This suggests that in such circumstances a realistic theory should lead 
to an h dependent result. 



9. Beyond LRT 

So far the elaborated spreading picture that we have introduced gave for D the same 
result as Kubo. So now we would like to see whether there are circumstances where 
this picture leads to novel physics beyond LRT. Considering s < and strong driving 
te < t^, it seems that one should substitute Ea. (fT4| in Eg.pSj). leading to a sub- linear 
dependence on the strength of the driving: 

D = h'^e^ t-\ [s<Ol[tp>t,] (18) 

However, one should be critical with regard to this result. The saturation of the 
coherent spreading process assumes a time- independent perturbation in Eq.®. This 
would be the case if cr = but not if ct = 2. 



10. Intrinsic dephasing hypothesis 

Considering the general problem of having a driving source with arbitrary spectral 
exponent cr, one realizes that as far as the integrand of the Kubo formula Eq.(IT]) is 
concerned the effective spectral exponent is s = sq — a. If a breakdown of this spectral 
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equivalence rule is observed, it constitutes a demonstration of a quantum anomaly. 
Indeed in the numerical experiment of Fig. 1, we contrast the dynamics which is 
generated by a time-dependent {a= — 2) perturbation, with that of a frozen (cr=0) 
perturbation that has the same s. Wc find that only the latter shows the coherent 
saturation of Eq. (fT4|) . 

Consequently we would like to conjecture that the implicit time-independence of 
the perturbation in Eq.Q leads to an intrinsic dephasing time which is finite in the 
limit — oo. We have in our problem only one time scale, the generalized Wigner 
time, and therefore the natural speculation would be ^ te- If this speculation is 
true, then the replacement in Eq. lfTS)) leads to the universal result Ea. ([TO)) . 

For the standard Ohmic case this implies that diffusive spreading persists beyond 
tg, and that LRT-like result oc still holds. For the non-Ohmics case this implies 
that the sub-diffusive or super-diffusive coherent spreading turns into normal diffusion 
with non-linear dependence on . Both expectations are supported by the numerical 
experiment of Fig. 1 that we further discuss and analyze below. 

11. RMT numerics 

In order to have a model that captures and tests the question of spectral equivalence 
it is convenient to use not the standard Wigner model, but rather its parametric 
invariant variation: 

U^{t) = + cos[/(t)] V'>^ + sin[/(t)] TAf (19) 

where V^^^ and V^^ are two independent realizations of a banded matrix, that has a 
bandprofile A|a;|^^~^, in the sense of Eq.(|3]), but with i. The energies En{t) are 
obtained via diagonalization of and the perturbation matrix Vnm{t) should be 
written in the same basis. The bandprofile of Kim is related to that of Vij as discussed 
in [11]. We have set g=l and A=l and verified that so ~ s\. We consider DC driving 
/ = const, and accordingly = A/^. 

In the time dependent adiabatic basis the perturbation matrix in the transformed 
Hamiltonian Eq.® is Wnm{t) = ifVnml{En—E„i), whose bandprofile is characterized 
by s = sq— 2. This matrix changes with time but it preserves its statistical properties. 
The question is whether its implicit time dependence generates an effective dephasing 
process. 

The numerical experiment is simple. On the one hand we make simulations with 
the time dependent Hamiltonian %. On the other hand we use a frozen version of 
Eq.® which we write in the ij basis as 

n,, (frozen) = £Aj + U,, + f W,, (0) (20) 

The initial state is assumed to be localized at i = 0, and an ensemble average over 
realizations is taken. Comparing the simulations (Fig.l) we deduce that there is 
intrinsic dephasing due to the implicit time dependence of the driving in Eq.®. 

In order to figure out what is the intrinsic dephasing time we plot in Fig. 2 
the scaled diffusion D/D^ versus the strength of the driving. The Ohmic case as 
conjectured is 'boring'. In contrast to that the sub-Ohmic and the super-Ohmic case 
exhibit departure from the universal expectation for large and small e respectively, 
indicating that the effective dephasing time becomes shorter than t^. We associate this 
systematic deviation with the infrared and ultraviolet cutoffs respectively: otherwise 
dimensional analysis implies that such deviation cannot emerge. We explain this 
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sensitivity as follows: The value of D is most sensitive to in the slow diffusion 
stage. In the sub (super) Ohmic case the slow diffusion stage is for long (short) times, 
and accordingly the sensitivity is to the lower (upper) cutoff, in-spite of the fact that 
the spreading is dominated by the high (low) frequency transitions. 

12. Conclusions 

LRT gives a finite classical-like result for the response of a low frequency driven Ohmic 
system. But in the case of a sub-Ohmic or super-Ohmic system, classical LRT predicts 
in the same limit either zero or infinite response. This is the notch where quantum- 
mechanics becomes most relevant, leading to an anomalous non-linear response. 

The analysis highlights the role which is played by the generalized Wigner time of 
|10| . which is the only relevant time scale in the universal continuum limit, and leads 
to a single-parameter expression Eg. dTUl) for the diffusion in energy space. 

The results might have a direct application concerning the heating rate of cold 
atoms in vibrating traps jl2j . where the experimentalist has control over both the 
shape (hence C{lu)) and the power spectrum (5(a;)) of the driving. In particular we 
note that if cold atoms are "shaken" without deforming the shape of the billiard, the 
spectral function that describes the fluctuations becomes sub-Ohmic fTSl. 
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Figure 1. (color online) Variance of energy spreading versus time for three 
different so-values, and two different driving rates /=5 (thick black lines) and 
/=12 (thin blue lines). The time-axis is rescaled according to the Wigner time 
tc given by Eq. l[7ll. Solid lines correspond to the simulations based on the 
Hamiltonian of Eg JlQII . while the dotted are for its frozen version Ea.l l20ll . 
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Figure 2. (color online) Dependence of diffusion on the rate of the driving for 
various values of so- The axes are X = e2/(2— s) ^^^^ y = D/e^^^'^~''\ where 
s = so— 2. Note that u)q < X < ujai where the level spacing is = 1 and the 
bandwidth is ui^i = 50. The deviation of D from universlity is due to the finite 
infrared or ultraviolet cutoffs: We see that in the superOhmic case (so > 1) the 
diffusion D becomes cjq independent for large e, while in the subOhmic case 
(so < 1) it becomes uici independent for small e. 



